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Abstract 


Designing an aircraft without conventional control 
surfaces is of interest to aerospace community. In 
this direction, smart actuator devices such as syn- 
thetic jets have been proposed to provide aircraft 
maneuverability instead of control surfaces. In this 
article, a numerical study is performed to inves- 
tigate the effects of unsteady suction and blow- 
ing on airfoils. The unsteady suction and blow- 
ing is introduced at the leading edge of the airfoil 
in the form of tangential jet. Numerical solutions 
are obtained using Reynolds- averaged viscous com- 
pressible Navier-Stokes equations. Unsteady suc- 
tion and blowing is investigated as a means of sep- 
aration control to obtain lift on airfoils. The ef- 
fect of blowing coefficients on lift and drag is in- 
vestigated. The numerical simulations are com- 
pared with experiments from the Tel- Aviv Univer- 
sity (TAU). These results indicate that unsteady 
suction and blowing can be used as a means of sep- 
aration control to generate lift on airfoils. 
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drag coefficient 
pressure coefficient 

fluxes of mass, momentum and energy 
viscous terms of the Navier-Stokes 
formulation 

non-dimensional actuator frequency, fc/U 0 0 

width of the orifice, m 

identity matrix 

free-stream Mach number 

transformation matrix from conserved 

variables to primitive variables dQ/dq 

conservation variables 

Prandtl number 

Reynolds number, /??9QooC 

A* oc 

actuator Strouhal number, fH /Vo 

temperature, °C 

contravariant velocities 

spatial mean jet velocity from actuator 

free-stream velocity, m/s 

speed of sound, m/s 

chord length, in 

total energy 

actuator forcing frequency, Hz 
pressure, nondirnensionalized by 7poo 
primitive variables 
distance outward from body 
time, nondirnensionalized by cja ^ 

Cartesian velocities in x and y direction, 

respectively, nondirnensionalized by a 

RMS velocity of the actuator jet oscillations 

Cartesian coordinates 

time step (nondimensional) 

angle of attack (degrees) 

difference operator 

ratio of specific heats 

coefficient of bulk viscosity 

viscosity 

density 

viscous stress tensor 
curvilinear coordinate directions 
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Subscripts 

max maximum 
min minimum 

x, y denotes differentiation in x and y directions, 
respectively 

ry, C denotes differentiation in // and 
C directions, respectively 
()oo free-stream condition 

Superscripts 

n denotes time level 

Q denotes quantities in generalized coordinates 

— denotes positive and negative flux 
conditions 

I. Introduction 


The experimental and computational investiga- 
tions of active control of flow past airfoils at high 
angles of attack is an area of active research as ex- 
tending the usable angles of attack has many im- 
portant applications. There are a number of arti- 
cles showing the effectiveness of flow control for air- 
foils. For example, in [10] leading edge suction was 
investigated for transition delay, in [23] jet flaps 
were employed for lift increase and in [6] surface 
suction/blowing was used to rapidly change lift and 
drag on rotary wing aircraft. However, most of the 
control techniques considered in the past required 
relatively more power input or involved weight and 
complexity penalties. 

In [1], [2] and [18], an innovative method for ac- 
tive control has been experimentally demonstrated 
using the so-called synthetic jet. The synthetic jet 
actuator produces a high-frequency jet from the 
surrounding fluid with zero net mass input. A 
novel feature of this actuator is that it requires 
only electrical power. In [1] and [18], a synthetic jet 
was used to produce a larger jet and in [2] a pair 
of actuators were used to show significant lift on 
cylinders. In another experimental work, Seifert 
et al [15] have investigated unsteady suction and 
blowing on a symmetric airfoil to increase post- 
stall lift. They show by introducing an unsteady 
jet near the leading edge tangential to the surface 
of a NACA0015 airfoil significant increase in lift 
can be obtained with relatively low momentum in- 
put. They also observed re-attachment of the flow 


and elimination of large wake region above the at- 
tached region. Later Seifert et. al. [16] and Seifert 
and Pack [17], investigated the dependence of ac- 
tuator location, momentum coefficient of the jet 
and frequency of the oscillation, and performance 
in higher Reynolds number. 

Some of the recent numerical simulations re- 
ported in [21, 6, 5] also support these findings. 
In [21], a periodic blowing and suction normal to 
the surface was used at 2.5% chord from the lead- 
ing edge of a NACA0012 airfoil. They showed 
that lift was increased for angles of attack be- 
tween 18°-35°. They used the Reynolds average 
Navier-Stokes (RANS) approach with the Baldwin- 
Lomax algebraic turbulence model. Hassan et. al. 
[6] have also used RANS approach with Baldwin- 
Lomax turbulence model. In their work zero net 
mass suction /blowing was placed at 13% chord and 
they found out that for certain oscillation frequency 
and peak amplitude, the lift can be increased albeit 
with high momentum input. The comparison with 
experiments were not given in both of these works. 
Finally, Donovan et. al. [5] reported numerical in- 
vestigations of both steady and unsteady jet on air- 
foils. They used an unsteady RANS incompressible 
flow solver with Spalart-Allmaras [19] turbulence 
model and compared their results with the experi- 
ments of [15]. Performance improvements were ob- 
tained by placing the actuator near the airfoil lead- 
ing edge. A significant lift increase of about 29% 
was obtained using synthetic jet actuators in the 
post-stall regime. However, their simulations as- 
sumed incompressibility which is not not strictly 
valid as the Mach number in the experiments of 
Seifert [15] was M =0.15. 

In this article, we present a numerical investiga- 
tion of unsteady suction and blowing on an airfoil 
using a compressible flow solver CFL3D [9]. We 
build upon the previous work of Donovan et al [5] 
to validate CFL3D for active flow control applica- 
tions and to address some of the anomalies in the 
previous reports. 

The plan of the paper is as follows. In §11, we will 
present the governing equations. In §111, we will 
present the computational methods used to solve 
them. In §IV, we will present numerical results for 
the baseline case. In §V, we will present numeri- 
cal results for the active control case. In §VI, we 
conclude the paper with a summary. 
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II. The Governing Equations 

The governing equations considered are the time- 
dependent, viscous compressible Navier-Stokes 
equations. The non-dimensional form of these 
equations in Cartesian coordinates for an ideal gas 
are 


dQ | dF j dG 1 f dF v , 
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where c p is the specific heat at constant pressure, 
and k is the coefficient of thermal conductivity. 
The pressure is defined by the equation of state 
for an ideal gas: 


P= (7 - l)[e- tj(u 2 + v 2 )}, 


(2) 


where 7 = and has a value of 1.4 for air. The 

<-v 

speed of sound, a, is defined as a 2 = 7^. The 

Reynolds number is defined as Re — y where 

c refers to a chord length and the subscript, 00, 
refers to free stream quantities. The dynamic vis- 
cosity, p, is approximated by Sutherland formula 


and equation (1) is closed by Stokes hypothesis for 
bulk viscosity (A 4- 2///3). 

A. The Coordinate Transformations 

In order to apply the numerical algorithm and 
boundary conditions easily, the governing equa- 
tions which are derived in the Cartesian coordi- 
nates, (x,y), must be transformed to the compu- 
tational domain or generalized coordinates, 

The transformation from Cartesian coordinates to 
general curvilinear coordinates for two-dimensions 
are 

T = ti £ = £(£, 2 /,*), r] = T){x,y,t). 

These coordinates conform to the surface of the 
body and maps the original Cartesian space (x,y) 
or physical domain onto a computational domain 
(£, //), which is rectangular with a uniform mesh. In 
terms of these curvilinear coordinates the Navier- 
Stokes equations become 
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where U and V are the contravariant velocities. In 
curvilinear coordinates the viscous flux terms are 
given by 

F v = J~ l (€xFv+£yG v ) and G v = J~ l (r} x F v + 7 f y G v ), 

where J — . Similarly, the stress terms in 

F v and G v are also transformed. 
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III. Computational Algorithm 


In this study the algorithm used to solve the 
Navier-Stokes equations is the CFL3D code re- 
ported in [9]. CFL3D solves the time-dependent 
conservation law form of the Reynolds-averaged 
Navier-Stokes equations. The spatial discretization 
involves a semi-discrete finite- volume approach. 
Upwind-biasing is used for the convective and pres- 
sure terms, while central differencing is used for 
the shear stress and heat transfer terms. Time 
advancement is implicit with the ability to solve 
steady or unsteady flows [13], [14], [20], [12]. Multi- 
grid and mesh sequencing are available for conver- 
gence acceleration. 


A. Time Differencing 


The Navier-Stokes equations (3) are discritized 
in time using the backward Euler implicit scheme 
and then the resulting nonlinear system is lin- 
earized in time about Q n to obtain 

[T + + S,,D n ]AQ n = R(Q"), (4) 
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Approximate factorization 


The central difference discretization to Equa- 
tion (3) results in a large banded square matrix 
which is sparse but would be computationally ex- 
pensive. To overcome this problem an approximate 
factorization is introduced which converts the two- 
dimensional operator into one-dimensional opera- 
tors: 
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This factorization yields tw T o block tridiagonal in- 
versions for each sweep. Prior to the execution 
of Equation (5), the corrections are constrained in 
order to maintain the positivity of the thermody- 
namic scalars p and p. For example, the update to 
pressure is taken as 


p n+i — p n + A[1 + + 


Ap 

pH 
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whenever < a c , wdiere a — —0.2 and <j) c — 2.0. 


B. Spatial Differencing 

The spatial differencing is done using a second- 
order-accurate upwind-biased scheme. The fluxes 
F, G representing pressure and convective terms 
are differenced using up-w T inding and a flux-vector- 
splitting method. For example, the flux difference 
in the £ direction is 

<%F = 3~F + +d+F-. 

The flux vector splitting is due to Van Leer. The 
flux F, for example, is split according to the con- 
travariant Mach number in the £ direction, see [9]. 
The split-flux differences are implemented as a flux 
balance across the ith cell holding spatial indices j 
and k constant as 

d-F++d+F- 
= lF+(Q-) + F-(Q + )}, +1/2 

-[F + (Q-) + F-(Q + )) i - 1/2 , 

wdiere F+(Q~) i+ 1 / 2 denotes a forward flux evalu- 
ated using the metric terms at the cell interface 
(i -f 1/2), and state variables are obtained by fully 
upwind second-order interpolation of cell-centered 
variables 


Q7+1/2 = (3/2)Qi - (l/ 2 )C?i-l> 
Qt + 1/2 = ( 3 / 2 )Qi+i - ( 1 / 2 )Q i+3 . 
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The diffusion terms are differenced using a second 
order central differencing: 

(d$Fv)i = (F v ) i+h - 


C. Convergence Acceleration 

A sequence of grids Go? G i, G 2 , , Gjv is 

defined, where Gjv denotes the finest grid and 
coarser grids are formed by successively deleting ev- 
ery other line in all coordinate directions. The fine 
grid serves to damp the high-frequency errors while 
the coarser grids damp the low-frequency errors. 
The coarse grids are solved with a forcing function 
on the right-hand side, arising from restricting the 
residual from the finer grids. The forcing function 
is the relative truncation error between the grids, 
such that the solution on the coarser meshes are 
driven by the fine grid residual. A fixed cycling 
strategy ( W-cycle) is used for the results presented. 
The solution were smoothed on each grid through 
five steps before switching to the next mesh. 

D. Turbulence Models 

In order to predict turbulent flows by solving the 
Reynolds averaged Navier-Stokes equations, clo- 
sure assumptions must be made about the turbu- 
lent stress and heat-flux quantities. For separated 
unsteady flow computations using Reynolds aver- 
aged Navier-Stokes equations, choosing the appro- 
priate turbulence models is not trivial. Among the 
available models, the Spalart-Allmaras (SA) model 
has been shown to be effective for a variety of flows 
including 2-D separated airfoil flows, see [5] and 
[8]. The Spalart-Allmaras model is a one equation 
model which solves a single transport equations for 
a modified turbulent viscosity. 

VI. Baseline Simulations 


The airfoil configuration used here is the one used 
at Tel-Aviv University (TAU) for low-speed wind 
tunnel test. This airfoil, we call it TAU0015, has 
a 0.4% chord notch at the leading edge and a 3% 
chord thick trailing edge, otherwise is a NACA0015 
airfoil. In order to lay a single-block structured grid 
around the airfoil, the original airfoil was modified 
to smooth out the square corners at the leading 


edge actuator. Figure 1 shows this modification in 
comparison with the original TAU0015 airfoil with 
the leading edge actuator. The baseline simulations 
given in this section uses no actuation. The follow- 
ing parameters were used in the simulations which 
are the same as in the experimental conditions: the 
Mach number of 0.15 and a chord Reynolds num- 
ber of 1.2 x 10 6 . The computations were performed 
using a 417X129 C-grid shown in Figure 2 with a 
minimum normal spacing of 0.0000015c. 

The boundary conditions applied are no-slip with 
no normal velocity at the body surface. Along the 
far-held upstream and the circumferential bound- 
aries, a quasi-one-dimensional characteristic analy- 
sis is used to determine the boundary data, assum- 
ing free-stream conditions exterior to the boundary. 
Along the downstream boundary, first-order ex- 
trapolation of the conserved variables is used. The 
outer boundary is sufficiently far away, r m ax = 
12c, from the airfoil, thereby minimizing the effects 
of the outer boundary on the flow over the airfoil. 

Figure 3 shows the enlarged view around leading- 
edge showing the grid near the actuator. Figure 4 
shows the computed lift coefficients versus angle 
of attack computed with experimental data. The 
results show very good agreement before stall, but 
deviate from the experimental data at stall, similar 
to the results of [5]. A drag polar is shown in Figure 
4 for the baseline case which also shows good agree- 
ment with experiment at lower angles of attack but 
deviate from the experimental values for post-stall 
angles of attack. Similar observations were made 
in [5] and [7] using the same airfoil and turbulence 
model but with different numerical methods and 
models. 

In [7] an unstructured grid was employed to 
study the effects of smoothing the square corners 
at the leading-edge actuator. Both the original 
TAU0015 and the modified TAU 001 5; see Fig- 
ure 1, were studied for the baseline case. They 
found remarkable agreement between modified and 
original TAU0015 airfoil computations but origi- 
nal TAU0015 model stalled before the modified 
TAU0015 model resulting in some disagreement in 
the post-stall region. They concluded that the dis- 
crepancies between computed and experimental re- 
sults was due to the way aerodynamic properties 
were computed and thus a proper interpretation 
of the experimental results is important. For ex- 
ample, experimental measurements used discrete 
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pressure taps and hence the lift coefficient consists 
of the pressure contribution only and leading edge 
actuator region was not included in the computa- 
tion. This was verified in [7] by recomputing the 
lift without including the actuator portion at the 
leading edge. 

We note here that we did not go back and recom- 
pute lift and drag; rather because the comparisons 
in Figures 4- 5 are sufficiently good that the remain- 
der of the paper will focus on this conditions. 

V. Control Simulations 


In this section, we present some numerical simu- 
lations of synthetic jet control and compare them 
with experiments. The experiments that we se- 
lected here for comparison used a leading-edge un- 
steady jet tangential to the surface. 

A. Actuator Boundary Conditions 

In all the calculations the synthetic jet actuator 
is modified using suction and blowing type bound- 
ary condition by prescribing velocity at the surface. 
The velocity boundary condition is given by 

u(x,0,t) = Asin(ujt)f(x), 

where the amplitude A = yj >, 
the prescribed frequency of oscillation uj = 
(F-rV27rr), the spatial distribution /(x) = 1, 
is the non-dimensional frequency and the os- 
cillatory momentum blowing coefficient < C it > is 
defined as 

< C M >= 2(i5T/c)(< u jet >/<7oo) 2 . 

The same C-grid used in the baseline simulations 
was used here. The jet at the slot was resolved 
using a fine grid consisting of twenty grid points. 

B. Numerical Results 

The parameters Mach number and Reynolds 
number were taken to be the same as in the base- 
line case. The experiments included various blow- 
ing coefficients < >, with F + = 0.58. To 

assess the effectiveness of synthetic jets, various 
blowing coefficients were tested at fixed post-stall 


angle-of-attacks a = 22°, a = 24° and a range of 

< > were used in the simulations. The turbu- 
lence model used was again Spalart-Allmaras with 
the time step At = 0.00005. In order to reduce the 
computational time required to converge to a so- 
lution with blowing, the baseline case solution was 
used as the initial flow conditions for the blowing 
computations. 

The variation of AC/ with blowing coefficients 

< Cy > for a = 22° is shown in Figure 6. The 

incremental lift increases as the momentum blow- 
ing coefficient increase for both the experimental 
and computational data in much better agreement 
than reported in [5]. The computational results 
over /under-predict lift. This may be due to in- 
ability of the turbulence model’s to predict sepa- 
rated flows, inconsistency in the actuator geome- 
try or the lack of grid resolution. Figure 7 shows 
the analogous variation in AC/ with blowing coef- 
ficients < C M > for a — 24°. The computed AC/ 
increase smoothly with increasing < > unlike 

the a — 22° case. For the uncontrolled flow, a large 
region of separated flow was seen with two coherent 
structures on the suction surface. The application 
of control makes the flow more attached which is 
consistent with the works of [4] and [5]. 

Figure 8 shows surface pressure distribution for 
the baseline and controlled cases. Figure 9 and 
Figure 10 show the lift coefficient time histories 
with synthetic jet control as a function of non- 
dimensional time (—?*) for a = 22° and a = 24°, 
respectively. Note that approximately 40 nondi- 
mentional times were necessary to obtain a statis- 
tically stationary solution. 

VII. Summary 

A computational investigation of tangential un- 
steady suction and blowing for separation control 
on an airfoil has been performed. The effects of 
zero net mass suction and blowing on lift increase 
were studied. The computed results were compared 
with available wind tunnel test results to determine 
the accuracy of the computational results. The 
numerical solutions were obtained by solving the 
Reynolds averaged Navier-Stokes equations. A grid 
resolution study was conducted using baseline (un- 
controlled) case to determine the appropriate grid 
density. The computed baseline results agreed rea- 
sonably with experimental results. For the active 
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control case, various blowing coefficients were in- 
vestigated at two angles of attack, a = 22° and 
q = 24°. The computed results were compared 
with the experiments from TAU. The results show 
reasonable agreement with the trends observed in 
the experiment. 

In general, the computations showed that the lift 
increased as the blowing coefficient increased. In 
all the control simulations, the grid resolution in 
general and time step-size in particular was found 
to be critical. 

Computational simulations shows that, tangen- 
tial unsteady suction and blowing on airfoil can be 
used as a means of separation control to generate 
lift. 
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FIG. 1 Leading edge regions of TAU0015 
and modified TAU0015 airfoils. 



FIG. 2 TAU0015 airfoil 417 x 129 C-grid, 
Ar m in = 0.0000015c, rmax = 12c. 
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